Lyapunov statistics and mixing rates for intermittent systems 
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We consider here a recent conjecture stating that correlation functions and tail probabilities of 
finite time Lyapunov exponents would have the same power law decay in weakly chaotic systems. 
We demonstrate that this conjecture fails for a generic class of maps of the Pomeau-Manneville 
type. We show further that, typically, the decay properties of such tail probabilities do not provide 
significant information on key aspects of weakly chaotic dynamics such as ergodicity and instability 
regimes. Our approaches are firmly based on rigorous results, particularly the Aaronson-Darling-Kac 
theorem, and are also confirmed by exhaustive numerical simulations. 
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I. INTRODUCTION 

It has been well known since the seminal works of 
Sinai, Ruelle, and Bowen (SRB) ^^ that the strongest 
chaotic systems (Smale's axiom A and Anosov systems) 
have SRB measures with exponentially decaying mixing 
rates (see also Ref. [4]). For these systems, the difference 
between temporal and spatial averages is statistically de- 
scribed by a Gaussian distribution (the central limit the- 
orem) and the convergence of both averages towards a 
unique value is assured. On the other hand, there is a 
wide range of systems where mixing rates and other re- 
lated correlation functions decay as power laws. Such 
class of dynamical systems, dubbed weakly chaotic in 
the Physics literature, typically exhibits weak statistical 
properties when compared with the chaotic ones. Exam- 
ples of weakly chaotic systems include maps with indif- 
ferent fixed points piSl]j billiards [9], and Hamiltonian 
systems with sticky islands in phase space [13], among 
others. These systems have in common an intermittent 
dynamical behavior, exhibiting a transition from regular 
to chaotic regimes which has attracted the attention of 
physicists and mathematicians in the last 20 years. We 
recall that the so-called mixing rate of a pair of phase 
space observable functions (j) and '0 for maps of the type 



Xt+l = f{xt) 
is defined as being the correlation function 



CM,4>) = 



0(x)^(/*(x)) dti{x) 
(f>{x) dfi{x) / tp{x) d^{x 



(1) 



(2) 



where ^.{x) is the invariant measure under the map /. A 
map is said to be mixing if Cf — > as i — >■ oo for any 
pair of phase space smooth observables (0, ■(/')• 



•Electronic address: 
t Electronic address: 
■••Electronic address: 



carlosamado@ig. com . br | 
asaa@lme . unlcamp . br | 



Since correlation functions as the mixing rate ([2]) might 
characterize the transition from strong to weak chaos, 
the attempt to relate them to more fundamental dynam- 
ical quantities is a goa l that looks, at first glance, really 
promising. In Ref. [11[, for instance, the finite time Lya- 
punov exponents 



Kt{x)^\Y.^n\r{f\x))\ 



(3) 



fc=0 



were considered for one-dimensional maps like ([T]) for 
which Aoo > 0, i.e., for strongly chaotic cases. Essen- 
tially, they show that if there exist two positive constants 
Aq and 7 > 1 such that 

/•CO 

Aft(Ao) - / ?7(At)rfAt < ait~\ (4) 

where ri[Kt) stands for the distribution of finite time Lya- 
punov exponents for the system in question, then we also 
have the following upper bound for the mixing rates 



Ct(0,V)<a2t"<^-'\ 



(5) 
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of any pair of Holder continuous observables (</>, V")- In 
other words, they have proved rigorously that if the finite 
time Lyapunov exponents tail probability is bounded by 
t~^ for large i, then correlations will be also bounded 
asymptotically by t~'''^~^\ These important results have 
inspired a recent work [l2| in which a related conjecture 
is made for weakly chaotic systems, i.e., irrespective of 
having koa > 0. The main results of |12| can be summa- 
rized as follows. 

1. They argue that scrutinizing the way in which 
1 — Mt decays for large t provides an "extremely 
efficient way" of studying quantitatively the decay 
of correlations. 

2. They conjecture, in view of recent results of [11| . 
that the estimates of @ and ^ are not optimal, 
and the decay properties of 1 — Mt and of Ct — 
J (j)d^ J tpdiJ, should be, in fact, both polynomial 
with the same exponent. 



3. They check numerically such a conjecture for an 
one-dimensional intermittent map with two indif- 
ferent fixed points of the Pomeau-Manneville type 
[13] , for which the polynomial decay rates of corre- 
lations are known exactly. Two other types of two- 
dimensional intermittent maps are also numerically 
considered to support the conjecture. 

Here, we show that this conjecture is false by pre- 
senting an explicit class of counter-examples. We will 
consider a general class of Pomeau-Manneville maps [1J| 
and show that their Lyapunov exponents tail probabil- 
ity (m decay faster than any power law, whereas their 
correlations ^ do exhibit a power law decay. Our ap- 
proaches are firmly based on rigorous results, particularly 
on the Aaronson-Darling-Kac theorem [lj|, and are also 
confirmed by exhaustive numerical simulations. For all 
maps in this class, correlation functions decay slowly than 
the Lyapunov exponents tail probability, suggesting that 
bounds of the type (j?]) and (O can be physically relevant 
also for weakly chaotic systems. 



to the spatial average / d dfj,. Then, for almost all initial 
conditions x G [0,1], the local expansion rate At{x) in 
@ converges to the unique positive Lyapunov exponent 
Aoo as t — > oo. On the other hand, if t is finite, At{x) 
assumes different values depending on the initial condi- 
tion X. The corresponding probability density function 
ri{At) = ??(A,i) is given by 



r/(A,t)= / S{At{x)-A)d^l{x). 



For large t, r]{A, t) takes the scaling form fji 

r/(A,t)-ry(Aoo,i)exp[-ifi(A)], 



(7) 



(8) 



where 51(A) > is a concave function with minimum 
at fi(Aoo) = 0. Then we have Q.{A) ~ ci(A - Aoo)^ 
and ri{Aoo,t) ~ (cit/7r)^/^, with ci > 0. Now, a simple 
calculation by using Laplace's method leads to 



(9) 



II. LYAPUNOV STATISTICS 

Our counter-examples consist in a general class of 
Pomeau-Manneville (PM) intermittent dynamical sys- 
tems of type ([H) with / : [0, 1] —> [0, 1], where 



/(x) r^ x{l + ax^-'^) 



(6) 



for a; — > 0, with a > and z > 1. The global form of / 
is irrelevant, provided it respects the axioms of an AFN- 
system [l^. For maps of the type ^, x = is an indif- 
ferent (neutral) fixed point, i.e., /(O) = and /'(O) = 1. 
Such systems are known to have power law invariant mea- 
sures near their indifferent fixed points. More specifically, 
we have d^{x) = ijj{x) dx, where aj(x) ^ bx^^'"^ near the 
fixed point x = 0, with a = {z — 1)~^ pjj. As conse- 
quence, such systems have diverging invariant measure 
near this points for z > 2. Moreover, finite invariant 
measure (1 < z < 2) implies ergodicity and the usual 
Lyapunov exponential instability, whereas the diverging 
case (z > 2) implies nonergodicity and subexponential 
instability. We will consider each of these cases sepa- 
rately and show that the conjecture proposed in [l2] fails 
for both. 



where fip — ^(Aq). The decaying properties of Eq. Q 
are definitively different than those ones predict by [13] . 
In fact, one has 



Mt 



1 exp(— iloi) 



(10) 



for Aq 7^ Aoo and Mt{Aoo) ^1/2. It is important to 
stress that there are many rigorous results in the liter- 
ature establishing polynomial bounds for the decay of 
correlations of Pomeau-Manneville maps in the regime 
1 < z < 2, see Q and references therein. Most notably, 
Sarig [6] and Gouezel Q have achieved optimal polyno- 
mial bounds for such correlations in this re gim e. There- 
fore, contrary to the conjecture proposed in [l2[, polyno- 
mial decay of correlations can occur simultaneously with 
exponential decay of Lyapunov tail probability distribu- 
tions. In fact, any conjecture stating that Mt should 
decay as a power law is generically violated for ergodic 
regimes. For the PM map with 1 < z < 2, this is indeed 
predicted by Theorem 2 (exponential level I result) of [8|]. 



B. Subexponential instability 



A. Exponential instability 

Let us first consider the statistics of finite time Lya- 
punov exponents ([3]) for randomly distributed initial con- 
ditions X S [0,1], in the case of finite invariant mea- 
sure cases (1 < z < 2). It is well known that ergodic- 
ity properties can determine c omp letely such statistics. 
For instance, BirkhofF theorem [T^ states that, in an er- 
godic regime, the time average of an arbitrary observ- 
able function t?, t~^ Sfc=o ^(/'^(^))' converges uniformly 



Let us consider now the cases for which the invari- 
ant measure diverges locally at the indifferent fixed point 
a; = 0, i.e., z > 2. For such cases, the system typically 
exhibits a nonergodic behavior and, hence, time averages 
do not converge to a unique constant value. Nevertheless, 
the Aaronson-Darhng-Kac (ADK) theorem [H, [H [13] 
ensures that a suitable time-weighted average does con- 
verge uniformly in distribution terms towards a Mittag- 
Leffler distribution of unit first moment. More specifi- 
cally, for a positive function ■(? and a random variable x 
with an absolutely continuous measure with respect to 



the Lebesgue measure on the interval [0,1], there is a 
(return) sequence {at} for which 



* fc=0 "^ 



Mfi 



(11) 



for i -^ oo, where ^q, is a non- negative Mittag-Leffler ran- 
dom variable of index a S (0, 1) and with unit expected 
value. The return sequence jotj for PM systems like 1^ 
and < a < 1 is given by [ll, [2l| 



11 /a\" sin(7ra) 

at~T-l — ) 1 , 

a \aJ ira 



(12) 



for i — > oo. What the ADK theorem is really pointing 
out here is to the explicit necessity of dealing with finite 
time subexpoential Lyapunov exponents 



t\^)-^Y.^^\r{f'im 



(13) 



fc=0 



instead of the usual exponents ([3]) for PM systems of the 
AFN type (see also Ref. ^2^]). From Eq. ([TT]), we have 






Saj 



(14) 



for t ^- oo, where the ADK average value (A) is given by 



1 / a \" sin(7ra) 



(A)^_(-j iriil^ / ln\f'ix)\ojix)dx. (15) 
ba \a/ na Jq 



The ADK theorem completely determines the correla- 
tions and the tail probability of Lyapunov exponents for 
the maps of the type ([5]), as one can see by considering 
a randomly distributed initial condition x G [0, 1] with 
probability density h{x) > in Eq. pTj) . leading to 



— ^f'd{f{x))h{x)dx= fi}dn 



(16) 



for i ^- oo. We can rearrange this expression and write 
Cticp,^)- f (t)dn f ddnr^a{d)t"'-\ (17) 



for t ^- oo, where (t>{x) — h{x)/uj{x) and the ADK aver- 
age {d) is given by an expression analogous to Eq. (1151) . 
It remains now to show that the tail probability of Lya- 
punov exponents for systems of the type ^ does not de- 
cay as predicted by Eq. ([T7|) . From Eq. P^ . one can ob- 
tain Mt for the map © by recalling that A^ = t^^"At, 
implying the following distribution of finite time Lya- 
punov exponents for systems of type ^ 



'7(At) = 



(A> 



■,('<■) 



(A) 



(18) 



where pa Ms a Mittag-Leffier probability density function 
with unit first moment, which corresponds to choice r" = 
aV{a), according to the definitions of [23|. Then, we have 
finally from Eqs. (g]) and (HH]) 



Mt 



it) 



p^:\s)ds, 



(19) 



for t — > oo, where u{t) = t^ "Ag/ (A). The behavior of 
Pa {x) for large x was recently discussed in [23], based on 
the known relation between Mittag-Leffler and one-sided 
Levy distributions 2J] and the Mikusinski's asymptotic 
analysis [25| of the latter. In particular, one has 



p'^:\^) 



A 2:(2"-l)/(2-2a) 



-cxp 



27ra 1 — a 

for r" = ar(a), valid for x — > oo, where 



(-Axi/(i^")), 



A- 



1-a 



a 



Via) 



i/("-i) 



(20) 



(21) 



The integral of Eq. (pn|) can be written in terms of the 
complementary error function, leading simply to 



Mt 
for large i, where 

B 



1 



erfc VBt 



2a 



1-a fAoTia) 



(A) 



i/("-i) 



(22) 



(23) 



The decaying properties of Eq. (P^ are also defini- 
tively different than those ones predict by Eq. (|17p . in 
the context of conjecture proposed in [12]. Once more we 
have 



Mt 



1 exp(-Bt) 



'27ra 



Bt 



(24) 



for large t, demonstrating finally that the conjecture pre- 
sented in [12| is false. It is noteworthy that Lyapunov 
exponents tail probability given by (pi| is essentially the 
same we would expect from finite measure cases, i.e., 
Eq. ((TU)) . This shows that the way in which Mt decays 
does not provide significant information on key aspects 
of weakly chaotic dynamics. 



III. NUMERICAL SIMULATIONS 

In order to test and illustrate the conclusions of the last 
Section, we perform an exhaustive numerical analysis of 
two particular AFN-maps of the type (]§]), namely the 
Thaler map [ia|, defined for z > 2 as 



/(■^) 



z-2 



„z-2 



-l/{z-2) 



(25) 



mod 1, and the modified Bernoulli map [22|, defined also 
for z > 2 as 

a; + 2^-la;^ < a; < -, 

f{x) = { 1 2 (26) 

x-2'-\l-xy, -<x<l. 



The Thaler map (|25p is very convenient here because 
its invariant measure density is explicitly known, namely 

M 



tu^x) 



.-1/" 



(1 



;-l/c 



(27) 



where a = (z— 1)^^, allowing in this way the explicit 
evaluation of the ADK averages like (|15p . In contrast 
with the Thaler map, there is no explicit expression for 
the invariant measure of the modified Bernoulli map, but 
it is known to have the form uj{x) ^ bk\x — Xfe|^^/", also 
with a = {z — 1)~^, in the neighborhood of each of the 
two indifferent fixed points xq = and xi = 1 [2g. Note 
that the ADK theorem is also valid for systems with more 
than one indifferent fixed point [l5|. For the Bernoulli 
map (pS)) . we also have return rates in the form at ^ i" 
for < a < 1 |15| . However, since its corresponding 
explicit expression for the invariant measure is lacking, 
we cannot evaluate the ADK averages for the modified 
Bernoulli map. We will show that this problem can be 
circumvented by exploiting the numerical data. 



A. Tail probability of Lyapunov exponents 

Our first task here is to determine if the tail probability 
of finite time Lyapunov exponents ^ for the maps ((25|) 
and (l^5|) do indeed decay as predicted by Eq. ([22) . From 
Sect. II, we known that such decaying behavior is assured 
if the the distribution of finite time Lyapunov exponents 
were effectively given by Eq. (fT8| . We compute numeri- 
cally the distribution of finite time Lyapunov exponents 
At for the maps ([25]) and (|26|) for random initial con- 
dition and large t, and confront the obtained numerical 
data with the theoretically predicted distribution (IT51) . 
The algorithm for the numerical computation of Mittag- 
Lefflcr distributions with arbitrary index a introduced 
in ^3] was instrumental to perform such task. The key 
point of our analysis is to check if a given distribution is 
well described or not by a generic Mittag-Leffler proba- 
bility density function. We recall that a Mittag-Leffler 
probability density pa (x) is defined from its Laplace 
transform as 



f 

Jo 



e-'"" p^;\x) dx 



^ r(i 

n=0 ^ 



{~sr°'y 



na) ' 



(28) 



for s > 0, with < a < 1. The choice r" = aT{a) 
assures that (x) =1, where the average here is evaluated 
with respect to pQ'(a;). From Eqs. (IT51) and (P5|) . we have 



A 

< 



< 

V 




FIG. 1: Graphics of (A^)/(A)^ as a function of a for Mittag- 
LefBer distributions with unit first moment, see Eq. (I30|) . The 
Mittag-Leffler index a can be determined from the value of 

(A2)/(A>2g[1,2]. 



the following constraints on the high order moments 



(A") _ nW'-^TiaY 
(A)" ~ nr{na) 



(29) 



of the probability density (fTS)) . One can evaluate (A") 
easily from the numerical data and the constraints (j29p 
can be objectively used to decide if a given distribution 
is well described or not by a Mittag-Leffler probability 
density. In particular, notice that one can determine the 
two free parameters of the distribution (fT8|) . (A) and a, 
by considering, for instance, (A) = i"^^(A) and 



(A2) aT{ay 



(A)2 r(2a) 



(30) 



It is very instructive to inspect the graphics of {A'^)/{A)'^ 
as a function of a, see Fig. [TJ For Mittag-Leffler dis- 
tributions with unit first moment, one has necessarily 
1 < (A2)/(A)2 < 2, with the boundaries correspond- 
ing, respectively, to a = 1 and a = 0. For such values 
of a, the Mittag-Leffler probability density function ap- 
proaches, respectively, a ^-function centered in a:: = 1 and 
a simple exponential e~^, see [2^. The violation of such 
boundaries would point out unequivocally that one is not 
leading with Mittag-Leffler distributions with first unit 
moment. Analogous bounds hold also for higher order 
moments ^, 1 < (A")/(A)" < n\. 

For the case of the Thaler map, both parameters (A) 
and a in the distribution (|18p are predicted theoretically 
by the ADK theorem, allowing the inspection of the con- 
vergence rate of Eq. (ITil) with respect to t and to the 
number of initial conditions used to evaluate the Lya- 
punov exponents. On the other hand, for the modified 
Bernoulli map one cannot determine exactly the average 
(A), but it is possible to infer its value by computing (A) 
from the numerical data and then using (A) = i^~" (A). 
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FIG. 2: Distribution of finite time Lyapunov exponents ([3]) for 
the Thaler map, determined from the iteration of Eq. (|25|l . 
with z = 22/7 [a = 7/15), for 2.5 x lO'^ initial conditions 
uniformly distributed on the interval [0,1]. The histograms 
are built directly from the numerical data, while the solid 
lines are the corresponding Mittag-Leffler probability density 
(|18|l , computed by means of the algorithm of Ref . |23l] . The 
inset and the background plots correspond, respectively, to 
t = 6 X 10* and f = 5 X 10^ 



Fig. [2] depicts the distribution of finite time Lyapunov ex- 
ponents for the Thaler map (PS)) . The plots show clearly 
that the distribution of Lyapunov exponents becomes 
peaked around the origin for i — >■ cxd, leading to (A") — >• 
for large t. In particular, one has Aoo = 0, in perfect 
agreement with the predicted distribution ((T5)) and the 
fact that the Thaler map is known to be weakly chaotic. 
Fig. [3] illustrates the case of the modified Bernoulli map 
(l26t . For all cases, we see, graphically and according to 
the higher order moments constraints ([2^ . that the dis- 
tribution of finite time Lyapunov exponents is very well 
described by a Mittag-Leffier probability density accord- 
ing to the prediction of Eq. p^ . The tail probability 
(E^ is then guaranteed for these maps. 

Our numerical examples are, in fact, illustrating the 
convergence of Eq. P^ , which is a consequence of ADK 
theorem. As expected, for large values of t and for large 
numbers of initial conditions, the histograms of both 
Figs. [5] and [3] approach the Mittag-Leffler probability 
density with the theoretical predicted values of a and 
(A). We could, however, detect another very interesting 
property. For a given value of i and a given number of ini- 
tial conditions, the corresponding histograms are already 
very well described by a Mittag-Leffler probability den- 
sity! With the increasing of t and the number of initial 
condition, such "instantaneous" Mittag-Leffler probabil- 
ity density approaches the ADK ones, as it is illustrated, 
for the Thaler map, in Fig. 2] and in Table ID The solid 
lines in both Figs. [2]and[3l for instance, are the instanta- 
neous Mittag-Leffler probability densities, i.e., their pa- 
rameters a and (A), although close to the theoretically 
predicted values, were calculated from the numerical data 



FIG. 3: Distribution of finite time Lyapunov exponents Q 
for the modified Bernoulli map, determined from the iteration 
of Eq. (ggi, with z = 28/13 (a = 13/15) and f = 6 x 10*, 
for 2.5 X 10^ initial conditions uniformly distributed on the 
interval [0, 1]. The histogram was built directly from the nu- 
merical data, the solid line is the Mittag-Leffler probability 
density computed with the algorithm of Ref. [23|. The be- 
havior of the distribution for large t is identical to the Thaler 
map case depicted in Fig. [3] In particular, one also has 
(A") — ^ for t — > CX3 and Aoo = 0, in agreement with the 
predictions of Eq. H18|) . 
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{A^)/(A)2 


1.305 
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1.300 
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(A'>/(A)' 


1.956 


1.947 


1.935 


1.921 


1.899 


(A*)/(A)* 


3.216 


3.189 


3.155 


3.117 


3.051 


(A^)/(A)5 


5.674 


5.598 


5.510 


5.414 


5.242 


(A«)/(A)« 


10.589 


10.383 


10.168 


9.927 


9.497 


15a 


10.81 


10.84 


10.88 


10.92 


11 



TABLE I: Statistical data for the graphics in Fig. [l] For 
each data set, the higher order moments constraints (|29p are 
respected, showing that each "instantaneous" histogram of 
Fig. Uis indeed well described by a Mittag-Leffler probability 
density. The values of the Mittag-Leffler index a for each one 
of the data sets (last row) were calculated from (A^)/(A)^, see 
Eq. (|30p and also Fig. [T] The last column corresponds to the 
probability density predicted from the ADK theorem, namely 
the curve (e) in Fig. [J] 

by using (A) = i^"" (A) and Eq. ^. Table U shows the 
values of the higher order moments constraints ([^^ for 
the data sets presented in Fig. |4l 



B. Correlation functions 

The numerical computation of the correlation func- 
tions ^ is rather tricky for the maps in question due 
to the highly discontinuous nature of the iterated maps 
p(x) for large t. For both cases (P5|) and (l^5|) . for in- 
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FIG. 4: Log-plot of the distribution of finite time subexpo- 
nential Lyapunov exponents (|13|l for Thaler map, determined 
from the iteration of Eq. ^, with z = 26/11 {a = 11/15), 
for 2.5 X 10^ initial conditions uniformly distributed on the 
interval [0, 1]. The set of points (a), (b), (c), and (d) corre- 
spond, respectively, to the histograms built from the numeri- 
cal data obtained for i = 10*, t = 5 x 10*, t = 25 x 10*, and 
t — 10®. Each one of these data sets are very well described 
by a Mittag-Leffler probability density, see Table |I1 The line 
(e) corresponds to the Mittag-LefHer probability density with 
the ADK values for a and (A). As one can see, the numeri- 
cally obtained distributions converges towards the prediction 
of the ADK theorem with the increasing of t. The increasing 
of the number of initial condition does not alter considerably 
such convergence, but a better description of the density tail 
requires a larger number of initial conditions, as expected. 



stance, the iterated map /*(a;) has 2* — 1 discontinuous 
points. An accurate numerical computation for large 
t of a correlation function as Ct would require an ex- 
tremely fine subdivision of the interval [0, 1], rendering 
the task practically and computationally unviable. Nev- 
ertheless, in order to establish the correlation function 
decaying (fTT)) . it is enough to assure that, for some value 
of < a < 1, the quantity 



t' = ft.Hn-)) 



(31) 



fc=0 



where i?(a;) is an integrable function, converges uniformly 
in distribution terms towards a random variable for large 
t. The ADK theorem does assure such a convergence 
with a — (z— 1)^^ for the maps in question, and, as 
demonstrated in Section 2, the correlation decaying (J17p 
is also firmly based on the ADK theorem. As an example, 
let us consider the observable function i9(a;) — sin^ ttx for 
the case of the modified Bernoulli map. According to the 
discussion of Section II, the correlation Ct{h,-d) for any 
smooth observable function h{x) > will exhibit a power 

lay decay as predicted by Eq. (fT7|) provided 91" does 
converge in distribution terms to a Mittag-Leffler random 
variable. Fig. [5] depicts such distribution and one can 
confirm the very good agreement with the predictions of 



FIG. 5: Distribution of the quantity 91 , given by Eq. (fSTjI . 
for &{x) = sin^ ttx, determined from the iteration of modified 
BernouUi map (gB)), with z = 28/13 {a = 13/15) and t = 
6 X 10*, for 2.5 x 10^ initial conditions uniformly distributed 
on the interval [0, 1] . The histogram is built directly from 
the numerical data and the solid line is the corresponding 
Mittag-LefHer probability density, computed by means of the 
numerical algorithm of Ref . [23] . 



the ADK theorem, assuring the validity of the correlation 
decaying (|17p. Similar results hold also for the Thaler 
map (P5|) and for other observables. 

We notice that the validity of Eq. (|T7)) is stronger than 
the ADK theorem, in the sense that the convergence to a 
Mittag-Leffler distribution given by Eq. (fTTT) is not a nec- 
essary condition to establish Eq. p^ . In fact, the exis- 
tence of a sequence at ~ t" such that a^^ X]fe=o '^if'^i^)) 
does converge in distribution terms towards a random 
variable, not necessarily of the Mittag-LefHer type, is 
enough to assure the decaying pT)). 



IV. FINAL REMARKS 

We close by noticing that the first map presented in 
[12| to support the conjecture we have just proved to be 
false is also a map with indifferent fixed points, namely 
the so-called Pikovsky map, which is defined implicitly 

by [27] 



/(.) + ^(l- 



0<.T<— , 

1 ^' 

f{x)Y, ^<2:<l. 



(32) 



The Pikovsky map is defined on the interval [—1, 1]. For 
negative x, one has simply f{x) = — /(— x). This map 
has two indifferent fixed points located at a; = ±1 for 
z > 1. The correlation functions for the Pikovsky map 
are known to decay as a power law [23, l28| . The authors 
of [12| present some numerical evidence suggesting that 
the tail probability Mt for the map (15^ would also decay 
with the same power law. This fact seems to contradict 



our results of Section II. However, a closer inspection 
of Eq. (15^ reveals that the Pikovsky map is not an 
AFN-map [15| and, hence, the ADK theorem cannot be 
invoked here to determine the distribution of finite time 
Lyapunov exponents. From the first equation of map 
([32|) . we have 



^ = (l-z)(2zx)-^/^ 



(33) 



showing that the axiom A (Adler's condition) [IQ] is not 
satisfied for a; = and positive z. The violation of Adler's 
condition here is related to the infinity slope of the map 
at the origin, and this is known to be capable of induc- 
ing some new dynamical properties as, for instance, the 



existence of a regular invariant measure in spite of the in- 
different fixed points, see Example 1 of [23, for instance. 
The failure of Adler's condition might explain why the 
authors of [IJl have arrived to the conclusion that Mj 
decays as a power law for the map (|32p. but certainly a 
deeper investigation of the Pikvosky map would be inter- 
esting and revealing. 
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